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Abstract 



Rotational friction on proteins and macromolecules is known to de- 
rive contributions from at least two distinct sources - hydrodynamic 
(due to viscosity) and dielectric friction (due to polar interactions). 
In the existing theoretical approaches, the effect of the latter is taken 
into account in an ad hoc manner, by increasing the size of the pro- 
tein with the addition of a hydration layer. Here we calculate the 
rotational dielectric friction on a protein {(df) by using a generalized 
arbitrary charge distribution model (where the charges are obtained 
from quantum chemical calculation) and the hydrodynamic friction 
with stick boundary condition, (C^^J^) by using the sophisticated the- 
oretical technique known as tri-axial ellipsoidal method, formulated 
by Harding [S. E. Harding, Comp. Biol. Med. 12, 75 (1982)]. The 
calculation of hydrodynamic friction is done with only the dry volume 
of the protein (no hydration layer). We find that the total friction 
obtained by summing up (df and Chyd^ gives reasonable agreement 
with the experimental results, i.e., (exp ~ (df + Chyd'^- 



*Email: bbagchi@sscu.iisc.ernet.in 



1 



1 Introduction 



In this article, we present an interesting result that the experimentally 
observed rotational correlation time of a large number of proteins can 
essentially be described as the combined effect of the rotational di- 
electric and hydrodynamic frictions on the proteins. Thus, one needs 
not assume the existence of a rigid hydration layer around the pro- 
tein, as is often assumed in the standard theoretical calculations of 
hydrodynamic friction. 

The study of rotational friction of proteins in aqueous solution has a 
long history [1 — 12]. Despite many decades of study, several aspects of 
the problem remain ill understood. For proteins and macromolecules, 
the rotational friction is obtained from Debye-Stokes-Einstein (DSE) 
relation given by, 

CR = 87Tr]R^ (1) 

where (r is the rotational friction on the protein and R is the radius of 
the protein. Naturally, the above relation assumes a spherical shape 
of the protein, which is often not correct. Moreover, there is ambigu- 
ity about the determination of some average radius of the protein. If 
one obtains the radius from the standard mass density of the protein 
(0.73 gm/cc), the values of the rotational friction are much smaller. 
The dielectric measurement of Grant |lj showed that the experimental 
value of rotational friction of myoglobin could only be explained by the 
above DSE equation, if one assumes a thick hydration layer around 
the protein, thereby increasing the radius of the protein. It is well 
known that spherical approximation embedded in DSE is grossly in 
error and the shape of the protein is quite important. However, even 
with the more recent sophisticated techniques such as tri-axial ellipsoid 
method \E\ and the microscopic bead modeling technique jEl El , which 
take due recognition of the non-spherical shape of the macromolecule, 
agreement with the experimental result is not possible without the 
incorporation of a rigid hydration layer [12]. In should be recognized 
that the effect of hydration layer thus introduced is purely ad hoc. 
In the case of tri-axial ellipsoidal method, the values of the axes are 
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increased proportionately by increasing the percentage of encapsula- 
tion of the protein atoms inside its equivalent ellipsoid fTU [T^ . On 
the other hand, the microscopic bead modeling technique uses beads 
of much bigger size P (3.0 A instead of 1.2 A) to take care of the 
effect of hydration layer. Without the hydration layer, the estimate of 
friction obtained from the theory is systematically lower. 

It has been recognized quite early that water in the hydration layer 
surrounding proteins and macromolecules has completely different dy- 
namical properties than those in the bulk The dynamics of wa- 
ter molecules in the hydration layer are also subject of great interest 
as they could play crucial role in the property and activity of these 
molecules. One often discusses the crossover from biological activity 
to the observed inactivity at low temperatures in terms of a protein- 
glass transition observed in the hydrated proteins fT^ . Recent inves- 
tigations have shown that the water molecules in the hydration layer 
are not only more structured but they also show slow translational 
and rotational motion than their bulk counterpart |T3, [THl HZI, [TSl HSI • 

Nevertheless, it is highly unlikely that water molecules in the sur- 
face of a protein such as myoglobin are so slow that we can replace it 
by a rigid hydration layer. On the contrary, all the recent experimen- 
tal and simulation studies have shown that the water in the surface 
of the protein exhibits bimodal dynamics j^. Majority of the wa- 
ter molecules seem to retain their bulk-like dynamics while a fraction 
(~ 20%) exhibits markedly slow dynamics. Recent solvation dynam- 
ics and photon echo peak shift experiment not only established the 
existence of slow water on the surface of proteins but also showed 
that the hydration layer is quite labile [21]. If one defines an average 
residence time to characterize the dynamics of water in the hydration 
layer, the residence time of bound or quasi-bound water is expected to 
range from 20 to 300 ps |22j . Question naturally arises how to under- 
stand quantitatively the role of the hydration layer in enhancing the 
rotational friction on the protein molecules. Clearly, the picturesque 
description of an immobile rigid layer around protein needs to be re- 
placed by a description where the hydration layer is slow but definitely 
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dynamic. 

This labile hydration layer has been explained in terms of a dynamic 
exchange model, which assumes that due to the presence of relatively 
stronger hydrogen bonding of water molecules with the charged groups 
at the surface of the protein, a surface water molecule can exist in 
either of the following two states - bound and free [231 • The free 
water molecules have dynamical characteristics similar to those of the 
bulk but the bound water molecules are essentially made static by the 
hydrogen bonding with the surface. In this picture, the slow time scale 
arises due to the dynamical exchange between the two states of the 
water molecules. Recent computer simulations seem to have confirmed 
the essential aspects of the dynamic exchange model (DEM) [24J . 

While the above model can provide a simple explanation of the 
origin of the observed slow dynamics, its correlation with dynamical 
properties of protein has not yet been established. This is a non-trivial 
problem as discussed below. 

The mode coupling theory (MCT) is another viable quantitative 
theory, which has been quite successful in describing translational and 
rational motion of small molecules [25j. This approach has also been 
extended to treat dynamics of polymer and biomolecules [2^. Let us 
recall a few of the lessons learned from the MCT of rotational friction 
of small molecules, and translational friction of ions in dipolar liquids. 
In both the cases, intermolecular dipole-dipole/ion-dipole correlations 
were found to play important role. It was also found that if one 
neglects the translational mode of the solvent molecules, then the fric- 
tion on polar solute increases by several factors. It should be noted 
here that the continuum models/hydrodynamic description of rota- 
tional friction always ignored this translational component. In fact, 
this translational component plays a hidden role in reducing the effect 
of the role of molecular level solute solvent and solvent- solvent pair 
(both isotropic and orientational) correlations that increase the value 
of the friction over the continuum model prediction. Thus, the issue 
is rather involved. In fact, the continuum model is found to give ac- 
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curate results due to cancellation of two errors: neglect of short-range 
correlations and neglect of translational contribution. In view of the 
above, it is thus important to note that the slow water molecules in 
the hydration layer can enhance the friction considerably. Thus, the 
classical picture of rigid, static hydration layer needs to be replaced by 
dynamic layer where the translational motion of the water molecules 
should be related to the residence time. However, only preliminary 
progress has been made in this direction. Thus, continuum models 
remain the only theoretical method to treat dielectric friction on com- 
plex molecules. 

An important issue in the calculation of the rotational friction is 
that proteins are characterized by complex charge distribution. The 
earliest models to estimate the enhanced friction on a probe, due 
to the interactions of its polar groups with the surrounding water 
molecules in an aqueous solution, employed a point dipole approxi- 
mation [1271 [2HI I29J . In the simplest version of the model, the probe 
molecule is replaced by a sphere with a point dipole at the center of the 
sphere. Such an approach is reasonable for small molecules, although 
continuum model itself may have certain limitations. The situation 
is quite different for large molecules like proteins because the charge 
here is distributed over a large volume and the surface charges are 
close to the water molecules. Thus, the point dipole approximation 
becomes inapplicable to such systems. This limitation of the early con- 
tinuum models was removed by Alavi and Waldeck |QOj who obtained 
an elegant expression for the dielectric friction on a molecule with ex- 
tended arbitrary charge distribution. By studying several well-known 
dye molecules, they demonstrated that the extended charge distribu- 
tion indeed has a strong effect on the dielectric friction on the probe 
molecules. The work of Alavi and Waldeck [20] constitutes an impor- 
tant advance in the study of dielectric friction. The role of dielectric 
friction has been studied for the organic molecules by other authors 

m 

The objective of the present work is to attempt to replace the rigid 
hydration layer used in hydrodynamic calculation. To this goal, we 
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calculate the hydrodynamic friction using the tri-axial method ||5j, in 
which the shape of a protein is mapped to an ellipsoid of three un- 
equal axes - closely representing the shape and size of the protein. No 
hydration layer is added in the calculation. We then calculate the di- 
electric friction using Alavi and Waldeck's model of generalized charge 
distribution for a large number of proteins. The friction contributions 
obtained from the above two methods are combined to obtain the total 
rotational friction. When compared, the total friction has been found 
to agree closely with the experimental result. 

We have also extended the work of Alavi and Waldeck to include 
multiple shells of water with different dielectric constants around a 
protein. The multiple shell model is introduced in concern with the 
experimental observation of varying dielectric constants of water from 
the hydration layer surrounding a protein to the bulk water. These 
shells have distinct dielectric properties - both static and dynamic. 
The resulting analytical expressions can be used to obtain quantitative 
prediction of the effects of a slow layer of water molecules on the 
dielectric friction on proteins. However, the multiple shell model in 
the continuum fails since it adds up the friction in every layer. This 
has been discussed in the appendix. 



2 Results and Discussion 

Below, we discuss the results obtained from the different aspects of 
rotational friction of proteins. The coordinates of the proteins are 
obtained from protein data bank (PDB) \^2\ ■ 

2.1 Dielectric Friction 

Dielectric friction is an important part of rotational friction for polar 
or charged molecules in polar solvent, because of the polarization of the 
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solvent medium. The solvent molecules, being polarized by the probe, 
create a reaction filed, which opposes the rotation of the probe. 

Many of the amino acid residues, which constitute the protein, are 
polar or hydrophilic. Therefore, in the aqueous solution, a protein and 
other polar molecules experience significant dielectric friction. There 
exist several theories (211 i28l [291 |33l El EH! , which account for the 
dielectric contribution to the friction. Some of these theories are con- 
tinuum model calculation of a point charge or point dipole rotating 
within the spherical cavity. Nee and Zwanzig provide an estimate 
of dielectric friction on a point dipole in terms of the dipole moment 
of the point dipole, dielectric constant of the solvent, Debye relax- 
ation time, and the chosen cavity radius. Later, Alavi and Waldeck 
\30\ extended this theory to incorporate the arbitrary multiple charge 
distribution of the probe molecule. 

The dielectric friction on the proteins has been calculated from the 
expression of Alavi and Waldeck for arbitrary multiple charge distri- 
bution model given below j2Dl, 

8 e, -1 - ^ /2/ + lx(/-m)! /ny/r,y 

Pricos Oi) Pricos 6j)cos{m(t)ji) 

where Rc is the cavity radius, (r^, (t)i) is the position vector and 
Qi is the partial charge of the ith atom. Pi^{cos{9i) is the Legendre 
polynomial. The maximum value of / used in the Legendre polynomial 
is 50. Cs is the static dielectric constant of the solvent. Since the 
solvent here is water, Cg is taken to be 78 and the Debye relaxation 
time td is taken as 8.3 picosecond (ps). 

The partial charges {qi) of the atoms constituting the proteins have 
been calculated using the extended Huckel model of the semi empirical 
calculation package of Hyperchem software. The dielectric friction is 
calculated on each of the atoms in a protein. The rotational frictions 
around X, Y and Z direction are calculated by changing the labels 
of the atom coordinates. The average dielectric constant C'jjp is the 



harmonic mean of the dielectric frictions along X, Y and Z direction. 
Here X, Y, and Z denote the space fixed Cartesian coordinate of the 
proteins, as obtained from PDB [32j. 

Table 1 shows the values of dielectric friction along X, Y, Z direction 
and their average. Continuum calculation method of the dielectric 
friction formulated by Alavi and Waldeck is dependent on the cavity 
radius and has been discussed in detail by them [|2Qj. They calculated 
the cavity radius from the observed orientational relaxation time of 
the organic molecules. The ratios of the longest bond vector of the 
organic molecules to the cavity radius ranged from 0.75 to 0.85. In 
Table 1, the calculation s of dielectric friction are performed using the 
cavity radius such that the ratio of the longest bond vector to the 
cavity radius is 0.75. 

In Table 2, we compared the average dielectric friction for the two 
above ratios - 0.75 (denoted as (^"p) and 0.85 (denoted as Cdf')- Cdf' 
is always larger than (jjp since the shorter cavity radius will put the 
charges close to the surface of the cavity, thereby increasing the polar- 
ization of the solvent and hence the rotational friction of the molecule. 



2.2 Hydrodynamic Friction 

The hydrodynamic rotational friction of the protein depends on its 
shape and size. Hydrodynamic friction was estimated earlier by the 
well-known DSE relation (Eq. 1). Perrin in 1936 ||36ij extended the 
DSE theory to calculate the hydrodynamic friction for molecules with 
prolate and oblate like shapes. Both prolate and oblate have two 
unequal axes. Harding |5j further extended the theory to calculate 
the hydrodynamic friction using a tri-axial ellipsoid. All the above 
theories employ stick binary condition to obtain the hydrodynamic 
friction. 

Tri-axial ellipsoidal technique requires the construction of an equiv- 
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alent ellipsoid of the protein. We have followed the method of Taylor 
et al. to construct an equivalent ellipsoid from the moment matrix 
P7j . The eigenvalues of this equivalent ellipsoid are proportional to 
the square of the axes. So this method provides with the two axial 
ratios. We then obtained the values of the axes using the formula 
given by Mittelbach ^ 



i?2 ^ (^2^52 + 0^) (3) 

' 5 

is the radius of gyration and A, B and C are the three unequal 
axes of a particular protein. 

Once the protein is represented as an ellipsoid with three principle 
axes, the hydrodynamic friction is calculated using Harding's method 
[SlEni- The hydrodynamic rotational friction of the ellipsoidal axes A, 
B and C are denoted as (a, Cb and Cc- The above rotational friction 
is obtained from the series of equations given below |3S] , 



Co = StttjABC, 

^""sABCiB^ai + C^as)' 

^^3ABC{A^ai + C^a:,)' 
2{A^ + B^) 



Cc = Ci 



ai 
ai 

A : 



3^50(^2^1 + 52^2) 

» dX 



(^2 + ^)^' 
dX 

(52 + A)A' 

00 dX 



(C2 + A)A' 
{A^ + X){B^ + X){C^ + X)f (4) 
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T] is the viscosity of the solvent. We have calculated the average of tri- 
axial hydrodynamic friction by taking a harmonic mean of the friction 
along three different axes, as given below, 



The values of hydrodynamic friction, along three principle axes (^4, 
B and C) of the ellipsoid and their mean, are tabulated in the Table 3. 
The A, B and C axes are not the same as the space fixed X, Y, Z Carte- 
sian reference frame. Note that the values obtained from the tri-axial 
method are much lower than the experimental values. Here, we can 
talk about an important aspect of standard hydrodynamic approach 
- hydration layer. One finds that hydrodynamic values of rotational 
friction underestimate the rotational friction unless the effect of hy- 
dration layer is taken into account. However, the effect of hydration 
layer is usually incorporated in an ad hoc manner, by increasing the 
percentage of encapsulation of the atoms inside the elhpsoid fEUHj- 
In this method, once the two axial ratios are obtained from the equiv- 
alent ellipsoid, the actual values of the axes are obtained by increasing 
the encapsulation of the protein atoms inside the ellipsoid. In the cal- 
culation presented here, the axes are obtained by equating with the 
radius of gyration. Therefore, we considered no hydration layer in 
this calculation of hydrodynamic friction. Later, we will show that 
this effect of hydration layer comes from the dielectric friction. 



2.3 Total rotational friction: Comparison with experimen- 
tal results 

We define the total rotational friction as the sum of dielectric friction 
{Cdf) Q-^d the hydrodynamic friction without the hydration layer (i.e. 
tri-axial friction, Cj.'^) as given below. 



1 Ir 1 1 1 1 



(5) 




/■ /-ati _|_ /-( 

'■^total — SDF ^ S[ 



TR 



av 



(6) 
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In Table 4, we have shown the values of the average dielectric {Cop), 
hydro dynamic {(tr) friction. Total friction {(total) defined above is 
shown in the fourth column. To compare with the experimental re- 
sults, we have shown the experimental values of the rotational friction 
in the next column. Note here, while the total friction, which is the 
contribution from both dielectric and hydrodynamic friction, is close 
to the experimental result, the microscopic bead modeling predicts 
the result, which is close to experimental value by itself The last 
column of Table 4 shows the references of the articles from which the 
experimental results are obtained. 

The similarity between the total friction and the experimental fric- 
tion is shown in figure 1, where we have plotted the experimental val- 
ues of rotational friction against the total friction for a large number 
of proteins. For most of the proteins, the results fall on the diagonal 
line. 

From the results shown in Table 4, we can conclude that the sum 
of dielectric friction and the hydrodynamic friction of the dry protein 
is approximately equal to the experimental results. 

Ctotal ~ Cexp (7) 



3 Conclusion 

Let us first summarize the main results of this work. We have cal- 
culated the hydrodynamic rotational friction on proteins using the 
tri-axial ellipsoid method, formulated by Harding [5], and the dielec- 
tric friction using the generalized charge distribution model derived 
by Alavi and Waldeck [30j. The hydrodynamic friction is calculated 
without the inclusion of any hydration layer. We have found that the 
combined effect of dielectric and hydrodynamic friction gives an esti- 
mate close to the experimental result. This approach seems to provide 
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a microscopic basis for the standard hydrodynamic approach, where 
a hydration layer is added to the protein in an ad hoc manner, to 
calculate rotational friction. 



The calculations adopted here are still not without limitations. The 
continuum calculation of dielectric friction is dependent on the as- 
sumed cavity radius. Unfortunately, there is yet no microscopic basis 
to assume certain value of the cavity radius for the calculation of di- 
electric friction. Moreover, the effect of increasing dielectric constant 
of the solvent from the vicinity of the protein to the bulk is not taken 
into account by Alavi and Waldeck j3Qj. Thus, we have attempted 
to incorporate a multi shell model to incorporate multiple shells with 
varying dielectric constants. The theory is described in the appendix 
in detail. The drawback of incorporation of multiple shells in the con- 
tinuum is that the frictional contributions from each of the shells add 
up, thereby giving rise to an unphysical large result. 



Similarly, the tri-axial method and bead modeling method suffer 
from the lack of microscopic basis to determine the exact values of the 
axes and the bead size, respectively. 



A potentially powerful approach to the problem is the mode cou- 
pling theory |10j , which uses the time correlation formalism to obtain 
the memory kernel of the rotational friction. The total torque is sep- 
arated into two parts - a short range part (which is called the bare 
friction Tfjare) and a long range dipolar part. The advantage of mode 
coupling theory is that it does not depend on any parameter. It uses a 
time dependent effective potential field in terms of density distribution 
and the direct correlation function given by [40j, 



Veff{r, t) = -keT J dr'dn'c{r - r', Q, n')Vnp{r\ n\ t) (8) 

The torque density is then expressed as, 

N,(r, t) = n(r, Q, t) [-VnVeff{r, Q, t)] (9) 

where, n(r,r2,t) is the number density of the tagged particle. The 
rotational friction comes from the torque-torque correlation function. 
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The final expression of ttie single particle (Fg) and collective friction 
(Fc) are given by [Hj, 



r oo r oo 



Tciz) = nare + A I e~^' I dk k' E F^k^t) cl^^{k) Fi,^{k,t) 



direct correlation function between any two dipolar molecules. Ff^^m 
and Fi^rn{k-,t) are the single particle and the collective orientational 
correlation functions, respectively. 

Eq. 10 and Eq. 11 are the standard mode coupling theory expres- 
sions for rotational friction. It has to be solved self consistently. In 
the overdamped limit, the self dynamic structure factor is expressed 
as. 



where fiim{k) = 1 — { — iy^{p/4:7T)ciim{k). I and M are the moment 
of inertia and the mass of the dipolar molecule, respectively. Tt{z) is 
the frequency dependent translational friction. 

The advantage of the mode coupling approach is that the once the 
charge density of the protein molecules and the dipole density of the 
water molecules surrounding the protein are defined, the rotational 
friction can be obtained in terms of the direct correlation function and 
the static and dynamic structure factors of the protein-water systems. 
These are again related by Ornstein-Zernike equation [42j. 




lil2m 



where, A = 
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The important aspect of this microscopic theory of dielectric fric- 
tion is the hidden contribution of the translational modes. In the 
hydration layer, the rotational friction is enhanced due to the slow 
translational component. This effect of translation could not be ap- 
proached through continuum calculation. Work in this direction is 
under progress. 



4 Appendix : Multiple shell model and the Draw- 
back 

Dielectric constant of water varies from the vicinity of the protein to 
the bulk water value. To understand the effect of this varying dielectric 
constant on the rotational dielectric friction of the protein, we have 
performed the continuum calculation of rotational dielectric friction 
using a multiple shell model. 

Nee and Zwanzig derived the dielectric friction contribution of a 
point dipole jSTj- Alavi et al. |3D] generalized it to obtain dielectric 
friction of a molecule with arbitrary distribution of charges. Castner 
et al. [43j generalized the point dipole approach to incorporate the 
discrete shell model with varying dielectric constant. Here, we have 
combined the approach of Castner et al. and Alavi et al. to obtain a 
generalized arbitrary charge distribution model for multiple hydration 
layers with varying dielectric constants around the protein. Figure 2 
shows the general scheme of this work. The protein is in the innermost 
cavity of radius a, where the water has a dielectric constant value of 4. 
The dielectric constant of water in the successive layers is assumed to 
increase up to the value of the bulk water, having a dielectric constant 
of 78. The width of each shell is assumed to be d. 

We first write down the electrostatic potential in two dimensions, 
which could be generalized to three dimensions using principle of su- 
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perposition. The electrostatic potential ^j{r) can be written as, 

$,(r) = $,.(r, 0) = f: Bl ^^^^ + A?, P,{cos 6) (14) 

where, j denotes the number of concentric shells surrounding the pro- 
tein. For n concentric shell, j can have a value from to n + 1. ji = 
denotes no boundary. The boundary conditions are, 

(i) Bf = Qivl, where Qi and are the partial charge and the position 
of the ith atom, respectively. 

(ii) $ ^ as r ^ oo 

(iii) $^(r^-) = ^j+i{rj+i), for j = 0, 1, 2...n. 

(iv) e,$;.(r,) = 6,-+i$;.+i(r,-), for j = 0, l,2..n, $;.(r) = 
(v ) A^' = 0, 

After incorporating the boundary conditions in Eq. 14, we get. 



i+i 



where, 



Bi = * ri {?i±lyTlU - ) (16) 



The reaction potential is given by. 



^jir, ^, 0) = E E Pi{cos-fi) (17) 

i=U=0 

where, 

Aqr m=l 

Plicosji) = —— E yr{o,, 00 (is) 

After few steps of algebra, we obtain the frequency dependent dielec- 
tric friction given below, 

N N 00 I 



\ 

\PJ^{cos9i)PJ^{cos9j) m cos{m<j)ji) x 



{I -my. 

{I + m)\ 
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E 

s=0 



Im 



a \2^+i / eg+i,4ma;) - 1 
efc,fc_i(ma;) + 



X 



(19) 



where, ejj_i = ej/ej_i, for all values of j. eg is the dielectric constant 
of the cavity. 



Above is the general expression of multiple (n) shell model. To 
write the final expression of dielectric friction for a two shell model, 
we assume Debye relaxation for the frequency dependent dielectric 
friction of two shells as given below. 



ei,o(ma;) = 1 + , (20) 
1 + I mu Tm 



e^Amuj) = 1 + ^ (21) 

1+2 mU TD2 

where, tdi and td2 are the Debye relaxation time for the first and 
second shell. 



The expression of dielectric friction for a two-shell model is given 
below, 



' 8^2/ + l^(/ 
Qdf = l^l^l^l^ 



m] 



tii=ii=im=ia W + 1 / (/ + m)! 
PricosOi) Pl^icosOj) cos{m(f)ji) x 



a 



X 



ei,o - 1 



a 



21+1 62 1 - 1 



+ 2 



a + dJ (262,1 + ly 

1 ^2.1 - 1 



a 



TD2 



a + dJ (261,0 + 1)2 (262,1 + ly 
(2ei,o + l)rD2 + (262,1 + l)rDi 



X 



(22) 
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The above expression has been numerically evaluated to find out 
the effect of dielectric friction on protein due to varying dielectric 
constant of water around the protein. The multiple shell model is 
found to overestimate the dielectric friction, as is evident from the 
above expression. 
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Table 1 



Table for the dielectric friction. The unit is 10^^ erg-sec. 
Cavity radius is chosen such a way that the ratio of longest bond 
vector {Rmax) of the protein to the chosen cavity radius {Rc) is 0.75. 



Molecule 


, — 

Rc (A) 




Cdf 


^ 

Cdf 


^ 

Cdf 


/-av 
^DF 


6pti 


29.50 


17.8 


13.2 


18.1 


16.0 


ligS 


26.10 


43.3 


36.6 


39.1 


39.5 


lubq 


34.30 


18.1 


18.3 


21.8 


19.3 


351c 


25.50 


52.3 


41.0 


41.9 


44.5 


Ipcs 


27.20 


90.5 


51.3 


66.1 


65.7 


lalx 


33.10 


63.0 


68.9 


49.5 


59.3 


Igou 


32.20 


43.8 


67.8 


103.6 


63.5 


laqp 


35.30 


44.5 


71.1 


132.1 


68.0 


le5y 


33.10 


98.9 


70.6 


89.9 


84.7 


Ibwi 


35.70 


78.3 


60.5 


108.1 


77.8 


IbSe 


33.50 


113.3 


112.2 


110.5 


112.0 


4ake 


50.30 


76.1 


170.8 


123.4 


110.7 


3rn3 


35.00 


118.8 


89.0 


56.8 


80.5 


Imbn 


28.00 


170.7 


162.0 


160.6 


164.3 
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Table 2 



Cavity size dependence of the dielectric friction. The unit is 
10~^^ erg-sec. 



Molecule 


/-0.75 


/-0.85 
^DF 


6pti 


16.4 


25.7 


ligS 


39.7 


61.3 


lubq 


19.4 


30.3 


351c 


45.1 


69.3 


Ipcs 


69.3 


111.0 


lalx 


60.5 


96.4 


Igou 


71.7 


114.6 


laqp 


82.6 


132.3 


le5y 


86.5 


136.4 


Ibwi 


82.3 


128.9 


IbSe 


112.0 


174.1 


4ake 


123.4 


211.7 


3rn3 


88.2 


138.1 


Imbn 


164.5 


263.1 


61yz 


107.8 


172.7 
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Table 3 



Table for the stick hydrodynamic friction using tri-8ixial el- 
lipsoid. The unit is 10^^ erg-sec. 



Molecule 


(A) 


— /'A 

Ctr 


/VD 

Ctr 


— /'C — 
Ctr 


/■av 
^TR 


6pti 


11.34 


57.8 


83.4 


85.1 


73.1 


ligS 


11.36 


72.9 


78.9 


84.9 


78.6 


lubq 


11.73 


71.2 


89.9 


94.0 


83.8 


351c 


11.51 


77.3 


84.5 


85.3 


82.2 


Ipcs 


12.38 


78.9 


106.5 


111.3 


96.6 


lalx 


13.47 


120.8 


127.3 


143.8 


129.9 


Igou 


13.61 


103.3 


141.7 


148.2 


127.7 


laqp 


14.45 


117.7 


171.1 


177.0 


150.1 


leSy 


13.81 


108.9 


145.7 


155.3 


133.4 


Ibwi 


13.94 


106.9 


155.4 


158.2 


135.7 


IbSe 


14.70 


167.5 


172.5 


178.2 


172.6 


4ake 


19.59 


298.3 


422.7 


442.8 


376.1 


3rn3 


14.31 


112.9 


166.5 


172.2 


145.1 


Imbn 


15.25 


163.7 


181.2 


210.1 


183.1 
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Table 4 



Comparison between the total friction and the experimental 
results. Results are given in the unit of 10~^^ erg-sec. The 

references to the experimental results of rotational diffusion of the 
corresponding proteins are given in the Ref. [8j . 



Protein 


rUB id 


^DF 


^TR 


Stotal 


sexp 


Bovine pancreatic trypsin inhibitor 


6pti 


1 n 
iO.U 


70 1 
(6.1 


on 1 


yo.o 


Calbindin D9k, holo form 


ligS 


oy.o 


( 0.0 


1 1 Q 1 

iio. i 


1 c; A 


Human ubiquitin 


lubq 


19.3 


83.8 


103.1 


118.9 


Ferricytochrome C551 


351c 


44.5 


82.2 


126.7 


130.1 


Plastocyanin, Cu(II) form 


Ipcs 


65.7 


96.6 


162.3 


149.5 


Oncogenic protein plS^^^*"^^ 


lalx 


59.3 


129.9 


189.2 


241.9 


Binase 


Igou 


63.5 


127.7 


191.2 


191.3 


Ribonuclease A 


laqp 


68.0 


150.1 


218.1 


186.1 


Azurin, Cu(I) form 


leSy 


84.7 


133.4 


218.1 


190.4 


Hen egg-white lysozyme 


Ibwi 


77.8 


135.7 


213.5 


203.6 


Bovine -lactoglobulin, monomer 


IbSe 


112.0 


172.6 


284.6 


270.6 


Adenylate kinase, apo form 


4ake 


110.7 


376.1 


486.8 


478.2 


Bovine Ribonuclease A 


3rn3 


80.5 


145.1 


225.6 


235.0 


Sperm Whale Myoglobin 


Imbn 


164.3 


183.1 


347.4 


246.3 
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Figure 1: The combined friction from hydro dynamic and dielectric is plotted against 
the experimental results. The sohd line shows the diagonal to guide the eye. 
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Figure 2: schematic diagram of the Molecular cavity and the hydration shell consti- 
tuted by the bound water molecules. The bulk water molecules are more randomly 
oriented 
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